AUTOMATED COUNTING of LEGO TOWERS 

Doron ZEILBERGER 1 

Dedicated to Henry Wadsworth Gould on his turning Fiq + Fj + F% years young 

Abstract: H.N.V. Temperley's method for counting vertically convex polyominoes is modified, 
generalized, and most importantly, programmed (in Maple). 

Preface 

I have never met Henry Gould in person, but have always admired his work. In particular, his 
charming article [G2] served me well many times, his (so far non-shaloshable) famous identities 
awed and frustrated me, and I still hope to use his work on Euler's constant[G3] to, who knows?, 
prove that it is irrational. I hope that this paper will please Henry, especially since his favorite 
sequence[Gl] makes a brief ('cameo role') appearance. 

Toys and Toy Models 

In spite of the many great triumphs of mathematics and science, there are many more problems 
that we can't solve than ones that we can. One of those, that so far defied us, is that of enumerating 
animals. Even Viennot's[V] powerful theory of heaps, that was so successful in enumerating directed 
animals (with the deceptively simple formula 3 n ), seems, at present, to be incapable of counting 
plain animals. 

A two-dimensional animal, alias polyomino, can be realized in terms of a LEGO tower. Suppose 
that we have an infinite supply of 1 x a (a > 1) LEGO pieces. Then every floor of the tower 
consists of a finite horizontal sequence of pieces separated by gaps. A vertical sequence of floors 
constitutes an animal if the resulting configuration is connected. Each floor can be described by 
the sequence of lengths of the pieces intertwined by the sequence of lengths of the gaps, so we 
have an infinite alphabet, each letter being of the form: a\, b\, 02, b<i, • • • > a™- Fixing the leftmost 
square of the bottom floor at the origin, one can view an animal as a word in this infinite alphabet, 
together with a specification of 'interfaces', which indicates where to place the leftmost square 
of the next floor in relation to the floor below it. The resulting creature is an animal iff the 
resulting configuration is connected. Note that in addition to the complication of having an infinite 
( in fact X^^Lo °° 2fc+1 ) alphabet, the condition of connectedness is global which makes it strongly 
non-Markovian, and hence so difficult. 

Whenever a problem seems impossible, it is not a bad idea to invent toy problems that are possible. 
Besides the positive chance that it might lead us eventually to solving the real thing, it is plain fun 
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to play, and to be able to get results. This is the case for polyominoes. 

The first toy model for animals was considered by Temperley ([T1],[T2] pp. 66-67) who treated 
vertically (equiv. horizontally) convex animals. These are LEGO towers in which every floor 
consists of a single piece. Since then, many other, much deeper, toy models were solved by the 
Ecole Bordelaise (e.g. [DV][V][B1][B2]), and the Australian school, under the doyenship of Tony 
Guttmann (e.g. [BGE]). 

In this paper Temperley's method is reviewed, its natural scope is realized, and then it is generalized 
in several directions. Everything has been programmed in Maple, and is contained in the package 
LEGO (available from my homepage), that can derive Temperley's generating function, and many 
others, instantaneously. 

Temperley's Method (slightly rephrased) 

Suppose that we have an infinite supply of 1 x a (a > 1) LEGO pieces. How many possible LEGO 
towers are there, with exactly one piece per floor, and with all the pieces parallel, whose total area 
is n unit squares? This problem is equivalent to that of counting vertically convex polyominoes, 
that was solved by Temperley in 1956. His method can be rephrased as follows. Let a{n) be the 
required number, then the generating function f(t) = Y^=i &{n)t n is the weight enumerator of all 
compositions (i.e. sequences of positive integers) with the weight 

r-1 

wt{a u . ..,a r ) = Y[{ai + Oi+i - 1) . 

i=i 

Indeed, every horizontally convex polyomino, of height r, gives rise to a composition (ai, . . . , a r ), 
where a, > 1 is the length of the i th floor, and to every pair of adjacent floors, of length a and b, 
there are a + b — 1 ways of placing the top one on top of the bottom one. 

More generally, let p(a, b) be an arbitrary polynomial, and L(a) an arbitrary affine-linear form 
L(a) = c\a + Co, where Co, c\ are integers, c\ > and c > 0. Define 

r-1 

wt(a u ...,a r ) = ^K)+-+£K) JJ p(o . ) a . +1 ) . 

i=l 

It is required to find 

/(t):=£>(C) . 

C*GC 

Here C is the set of all compositions. 

Temperley's trick, that is hereby promoted to method, is to consider the weight enumerators F(a) 
for the subset of C whose first component is a, 

F(a) := wt(a,a 2 , ■ ■ ■ ,a r ) , 

(a,02,...,a,.)€C 
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and consider the two-variable generating function 

oo 

^(z,t):=Y,F(a)z a . 

a=l 

Once we know <3? we would also know f(t), since f(t) = $>(1, t). 
The natural equations for the F(a) are 



F(a) = t L{a) + t L ^ J2 P(°> b ) F ( b ) > (*) 
6=1 

since every composition C, that starts with a, is either (a), whose weight is t L ^ a \ or is of the form 
C = (a, C), where C is a composition on its own right, whose first component starts with, say, b, 
for some b > 1, and then wi(C) = t L ^p{a,b)wt(C). 

Now let's expand the polynomial p(a, b) in powers of b: 

R 



P( a , b ) = ^2Pr(a)b r , 



r=0 



and plug it in (*), to get 

oo / R \ i? /oo \ 

6=1 \r=0 / r=0 \6=1 / 

Now multiply both sides by z a , and sum over a > 1, to get 

oo oo R / oo \ / oo \ 

J>(a)z a = ^^W^ P>> r (a)t L ^ HT W . (**) 

a=l a=l r=0 \a=l / \6=1 / 

Let's define, 

oo oo 

h(z,t) :=jN L(o) z a , and g r (t, z) := ^ Pr {a)t L ^ z a , (0 < r < R), 

a=l a=l 

which are certain explicitly computable rational functions in (t, z). Also define 
Eq. (**) now becomes: 

R 

Q(z,t)=h(z,t) + Y,g r (z,t)& r \l,t) . (***) 



r=0 
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Now apply (z-^) s to both sides of (* * *), for s = 0, 1, . . . , R, and then plug in z = 1, in order to 
get R + 1 linear equations, with coefficients that are rational functions of t, for the R + 1 unknowns 
<&[ r l(l,£), < r < i?, solve them, and get in particular $(l,i) = /(i). □ 

The procedure that implements this in LEGO is LEGO. After you downloaded LEGO to your working 
directory, get into maple and type read LEGO:. The function call is LEG0(L,p,a,b,t). For ex- 
ample, to get Temperley's original generating function, type LEG0(a,a+b-l,a,b,t) ;, and Maple 
would respond with 

*(* ~ I) 3 
At 3 - 7t 2 + 5t-l ' 

Another simple example (attributed to Moser in [K]) is to the enumeration of leftist (equivalently 
rightist) horizontally convex polyominoes. Here the leftmost point of every floor is farther to the 
left (or right on top) than the leftmost point of the floor below it. Here p = b and the function call 
is: LEG0(a,b,a,b,t) ;, producing the output: ^3^2 , which equals (why?) Y^=i F2n-it n , and 
hence the number of leftist horizontally convex polyominoes with area n is F 2n -\. 

To get the number of LEGO towers, with one piece per floor, and where every floor is perpendicular 
to the floor below it (so we have a kind of zig-zag pattern, that lives in three dimensions), do 
LEG0(a,a*b,a,b,t) ; , and you would get: 

t(l - 3t + It 2 - t 3 ) 
1 - 5t + 6t 2 - 3t 3 + t 4 

The reader is welcome to think up other kinds of towers, and to use LEGO to enumerate them. 
Building With Colored Pieces 

If we have s different colors of LEGO pieces, i = 1, . . . , s, and that the area of a piece of color i 
and length a is given by the affine- linear function Lj(a). Suppose also that the number of ways of 
placing a piece of color j and length b on top of a piece of color i and length a is Pij(a, b), where 
the pij(a,b), (1 < i,j < s) are polynomials in (a, b). 

This amounts to weighted counting of colored compositions. A colored composition is a sequence 
of colored integers (a[ kl \ a^ 2 " 1 , . . . , a\ ki \ ■ ■ ■ , a^^), where the superscripts denote the colors, (so 
1 < ki < s), and 

«*(a< fcl \a< fa \...,aW ...^ . 
It is required to find 

f{t) := £ wt{C) . 

C€CC S 

Here CC S is the set of all colored compositions with s colors. 
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For each color i (i = 1, . . . , s), and for each integer a > 1, consider the weight enumerators Fi(a) 
for the subset of CC S whose first component is , that is: 

F,(a):= Yl wt{a^\a { ^\...,a^)) , 

(a('),4 fc2) ,..,o^>)eCC s 

and define the two-variable generating functions 

oo 

$i(z,t) := ^2F l (a)z a (1 < i < s) . 

a=l 

Once we know the $i we would also know f(t), since f{t) = Yli=i 
The natural equations for the Fj(a) are 

oo s 

F l (a)=t L ^+t L ^Y,Y,P^ a ^ F ^ > (*) 

6=1 j=l 

since every colored composition in CC S , that starts with o^, is either (a^), whose weight is t Li< - a \ or 
is of the form C = (a' 1 ', C"), where C" is a composition on its own right, whose first component starts 
with, say, b^\ for some b > 1, and some color j (1 < j < s), and then wt(C) = t Li( - a ^pij(a, b)wt(C). 

Expand the polynomials Pi t j(a, b) in powers of b: 

Pi,j(a, b ) = J2ptj(a)b r , 

r=0 



and plug it in (*), to get 



J.r)fs h r\ _ + L t (a) , ,L 4 (a) V^JO, 



j = l 6=1 yr=0 y j = l r=0 \6=1 / 

Now multiply both sides by z a , and sum over a > 1, to get 

OO OO S / OO \ / OO \ 

E^(^° = E* Li(0) ^ + EE £p$(«)* Li(0) *° E^ r • (**) 

a=l a=l j=l r=0 \a=l / \6=1 / 

Let's define, 

oo oo 

h i (z,t):=J2t L * (a) z a , and <$(t,z) := ^ P Q(a)t L ^z a , (0 < r < i?^, 1 < i,j < s), 

a=l a=l 

which are certain explicitly computable rational functions in (t, z). Eq. (**) now becomes: 

S 

$ i (z,t) = h i (z,t)+Y,£&J(z,t)*j ] (l,t) ,(!<*<*)■ (***) 

j = l r=0 



Now apply (z-^) 1 to both sides of (* * *), for I = 0, 1, . . . , max.,- Rij, and then plug in z = 1, in 
order to get Xli=i( max i + 1) li near equations, with coefficients that are rational functions of 
t, for the unknowns $j i] (l,t), < I < maxj Rij, solve them, and get in particular $j(l,t), for 
i = 1, . . . , s, and finally /(t) = £- =1 $,(1, i). 

The procedure in LEGO that implements the weighted enumeration of colored compositions is 
muLEGO, the function call is muLEG0(Ls,p,a,b,t) ; , where Ls is the list [Li, . . . , L s ], and p is 
the list of lists 

[[Pl,l,Pl,2, • • • ,Pl, s ], • • • , \Pi,l, ■ ■ ■ ,Pi,j, ■ ■ ■ ,Pi, a ], [p s ,l, • • • ,P„,j, ■ ■ ■ ,Ps,s]\ , 

a, 6 are the variable names, and t is the variable chosen for the generating function. For example, 
to find the generating function for locally stable horizontally convex polyominoes, (i.e. the center 
of gravity of every floor is in the interior of the floor below it) do: 

muLEG0([2*a,2*a-l] , [[2*a-l,2*a] , [2*a-2 , 2*a-l] ] ,a,b,t) ;, getting the output t<yl ^~!i i ~f ) ■ 

Yet another example is the number of LEGO towers, with one piece per floor, but now you have 
an infinite supply of both 1 x a and 2 x a pieces. Furthermore, the towers are to be constructed in 
such a way that all the pieces are parallel to each other (they each have a designated length-side 
and width-side, even 1 x 1, 1 x 2, 2 x 1, and 2x2 pieces). The function call is 
muLEGOC [a,2*a] , [[a+b-l,2*(a+b-l)] , [2*(a+b-l) ,3*(a+b-l)]] ,a,b,t) ;. I omit the output 
(do it yourself!). More generally, to find the number of such towers where you have an unlimited 
supply of pieces of the shape: 1 x a, 2 x a, . . . , R x a, (where a > 1) , you may use the built-in function 
MIGDAL(R.t) ;. Thus Temperley's original generating function is, in particular, MIGDAL(l.t). 

Higher Dimensional Structures 

Suppose that we have an infinite supply of a x b pieces 1 < a,b, how many towers can we build 
of side-surface-area n (we don't count the area of the base and top) where each floor has exactly 
one piece, and all the 'lengths' are parallel to each other? Now we have vector compositions 
{[a^\ a^], [a^jOg 2 ^], • • • , [a^, Or 2 " 1 ]), where the weight is 

Wt ([a?\ a?'], [a«, af >], . . . , [a^,a^]) = ^+^) JJ^D + a V _ 1)(a? ) + a g) _ 1} . 

More generally, fixing m, and using vector notation a = (a^, . . . , c^™)), we have to weight- 
enumerate sequences (ai, . . . , a r ) with the weight given by 

r-l 

wt( ai , ...,a r )= ^K)+-+^K) Yl p(a t ,a i+1 ) , 

i=i 

where L(a) = L(a ( - 1 \ . . . , a^) is affine-linear in its variables and p(a,b) is a polynomial of 2m 
variables. The previous analysis goes almost verbatim, and is left to the reader. 
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The procedure in LEGO that handles this case is LEGOmul. The function call is LEGOmul (L , p , a , b , t) ; . 
Here L is the affine-linear form in the m variables a; p is a polynomial in the 2m variables a,b; a 
and b are the two lists of m variables, used to describe L and p; and t is the designated variable of 
the generating function. For example to solve the problem described above, type: 
LEGOmul (2*al+2*a2 , (al+bl-1) * (a2+b2-l) , [al,a2] , [bl,b2] ,t) ;. 

Building With Colored Multi-dimensional Pieces 

Suppose that we have s different colors, i = 1, . . . , s, where the pieces of color i are c^-dimensional 
pieces, ai xa 2 x... xa^, with the a\, ai, ■ ■ ■ , 0,4. > 1- The discussion on colored one-dimensional 
pieces goes almost verbatim, only now the a^, bi, and Zj are multi- variables (with di variables). 

The function call is: muLEGOmul (Dims, Ls.pols, a, b,t), where Dims is the list of dimensions of the 
colors i = 1, . . . , s, Ls is the list of affine-linear functions, where Ls [i] depends on a[l] , ... , 
a [Dims [i] ] ; pols is the list of lists of polynomials pij(di, bj), where a>i stands for (a[l], . . . , a[Dims[i]]), 
and bj stands for (6[1], . . . , b[Dims[j]]); a and b are the letters chosen to express the indexed vari- 
able; and t is the designated variable name for the argument of the output. 

Here are two examples: 

muLEGOmul ( [2 , 2] , [a[l]+a[2] ,2*a[l]+2*a[2]] , [ [(a[l] +b [1] -1) * (a [2] +b [2] -1) , (a[l]+b[lj- 
l)*(a[2]+b[2]-l)], [(a[l]+b[l]-l)*(a[2]+b[2]-l) , (a [1] +b [1] -1) * (a [2] +b [2] -1)] ] ,a,b,t); 

muLEGOmul ([1,2] , [a[l]+l,a[l]+a[2]] , [ [a[l] +b [1] -1 , (a[l] +b [1] -1) *b [2] ] , [a [2] *(a [1] +b [1] - 
l),(a[2]+b[2]-l)* (a[l]+b[l]-l)]] ,a,b,t);. 

Future Directions 

The present generalizations of Temperley's method should be extendible much further, for example 
to the enumeration of convex polyominoes, both according to area and perimeter (see [B1][B2] and 
references thereof). But, now we no longer get rational functions, and the natural context would 
be functional and functional-differential equations, that might, if in luck (like in [DV]) turn out to 
be an algebraic generating function, in which case the future program should be able to guess it 
empirically, and then prove it rigorously. 

Now the natural equations for the F(a) would be: 



b=l b=a+l 

where p(a, b) and q(a, b) are different polynomials of (a, b) or (q®, q% , . . . , q\, q%, ■ ■ ■)■ 
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